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545. Recurrente betrekkingen. 


(45.1) Bij de numerieke behandeling van differentiaalvergelijkingen 
loopt men voortdurend tegen recurrente betrekkingen aan. Daarom 
bespreken we nu eerst enkele eigenschappen ervan (een deel hiervan 
vindt men al in TA (12.9) t/m (12.16)). 


(u5.2) Onder een lineaire k+1=terms recurrente betrekking voor een rij 
tai);ao verstaat men een relatie 

(45.3) aisk + C4iâiak-1 + . . e + Criai = Biak? ll = EE 

met gegeven esi J = 1sveesks i = 0,1,2s... Men noemt zo'n rij 

{a;} ook wel een oplossing van (45.3). Een oplossing is blijkbaar 


bepaald door (willekeurig voor te schrijven) beginwaarden 


Agsajs: . arj 


(U5.u) ‘Een bijzonder geval is dat de esi niet van i afhangen, dus b.v. 
in welk geval men spreekt van een reeurrente betrekking met constante 
coëfficienten. 

(u5.,5) Als 8j = 0 voor alle j noemt men de recurrente betrekking 


homogeen. 


Men verifieert onmiddellijk: 
(u5.,6) Als (85.3) homogeen is vormen de oplossingen een k-dimensionale 
lineaire ruimte d.w.z. dat als (a, ji} en fa, j} twee rijen zijn die aan 
(65.3) met alle 8j=0 voldoen, dun Gout Sk elke lineaire combinatie 
{pa,j*aa,;)}; en er zijn rijen agih, -lajs} zodat elke oplossing 
van de homogene vergelijking een lineaire combinatie hiervan is 
(neem nl. maar voor fa, ;} de rij waarvan de eerste k termen luiden 
1,0,0,...,0, voor las} de rij die begint met 0,1,0,...,0 etc.). 
(H5.7) Een willekeurige oplossing van de inhomogene vergelijking 
(H5.3) wordt verkregen door een vaste oplossing van (45,3) te 
vermeerderen met een geschikte oplossing van de homogene recursie 
die bij (U5.3) past. Blijkbaar vormen alle oplossingen van (u5,3) 
een k=-dimensionale lineaire variëteit. 
(u5.8) Van een lineaire homogene recursie met constante coëfficienten 
laten de oplossingen zich op zeer eenvoudige wijze voorstellen: 


Stelling. Zij 


(u5.9) t Ciajikeg Pee + eKai = 0 


Liek 
een recurrente betrekking met constante coëfficienten. 
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Als de „eakteristieke vergelijking van deze recursie: 


(45.10) (@Wexk + og 1 ve te = 0 


uitsluitend enkelvoudige wortels Ar see e sg heeft dan vormen de 
rijen ut 5) i>0o een basis. van de oplossingsruimte, d.w.z. voor een 
Wil ekenrizs am {a; } bestaan constanten Qi s--.>d zodat 
di 
(u5.11) aj s Gale WM want Aer 
Als, algemener, de karakteristieke vergelijking de verschillende 


wortels Ar eee) heeft met multipliciteiten Mi sees dan vormen 
de rijen CiTG), i>g, Mm = 0, shar eee ie j = 1,...sps een basis van 
de oplossingsruimte, dew.z. 

(45,12) a; = qilidAi +. + ap (aj 


met a; een veelterm van ED en 
Bewijs. Zij im 5 i(i-1)...(i-mt1) voor m geheel > 0, io= 1, voor 
alle i. Wegens ema, ) = 0 voor m = 0,1,...sl1=1 geldt 


mm 8 
TROEL en, = 0 (regel van Leibniz voor het herhaald differen- 
tieren van een produkt). Dus 

[Ciek) en + ec, Citk- 1) Pii tees + Cr ÂnXx ij 5 fs 


Hieruit wat meteen dat {(i Ai Ra i>oe f 5 0O,1s...shi=1 een oplossing 
is van (45,9). Hetzelfde ale voor À3 seenshge 

We krijgen op deze manier precies k oplossingen, en zij zullen 
dus een basis van de oplossingsruimte vormen als ze lineair onafhan- 
kelijk zijn. Dat dit laatste inderdaad het geval is laten we een-- 
voudigheidshalve aan een speciaal geval zien; het algemene geval 
gaat echter analoog. Zij p=3, Ui1=1, 4222, U3=3. We noteren van 


elk der genoemde oplossingen de eerste 6 termen. 


1 A, Ai M Mi 
1 A MM Mî x$ 
0-1 2À, 3MÂ uag 5M% 
1 As MA A5 LE; A5 
0 1 2As 3M3 A3 5A$ 


0. 0 2.1 3.23 U. 5.uM 

We beweren dat reeds deze 6 rijtjes lin-onafhankelijk zijn en 
daarmee de hele oplossingen. De onafhankelijkheid van deze rijtjes 
volgt uit die van de kolommen. Stel derhalve dat ri: Xx le kolom + 

+ ra Xx 2e kolom + ... = 0. Zij g(x) = Pit Pix +... + re x° . Dan 
geldt blijkbaar g(A,) = g (A2) = gl) = glAs) = g'(A2) = B"(A3) = 
0, zodat g deelbaar is door Ox-A, Gerda DP Ox-A, WS en dat kan wegens 
de graad van g alleen als g identiek 0 is, dus nz As ee == Fez 0. 
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Rest nog te bewijzen dat ook Et Fe i>o, M= 0, speen oje1 
j = 1s...,p een basis vormen, Dit VoleE uit het feit dat ook dit 
k oplossingen zijn en dat de eerder genoemde basis hieruit door 
lineaire combinatie ontstaat, zulks omdat {li} lin. combinatie 


is van (if), [im1j, amen he 


(u5.13) Men kan een (k+1)-terms recurrente betrekking vertalen in een 


(U5,15) 


twee terms vector recurrente betrekking, hetgeen voor het vervolg 


heel nuttig zal blijken. Definieer nl. rijen vectoren {vv} en {b;} 


en matrices A; als volgt: 0 4 En 
0 0 0 Ae eer es 
Te , OE: 
(u5.14) Vs = 5 n b. = kg Â. z= d Ln , d 
ì k L 6 ì 6 0 n el 
itk-1 ne zie ER 
Bi+k ki k-1,i 1,i 


Dan wordt (45,3) 


Baar Av; + b;» GEK 0 Dr 


(U5.,16) In het bijzondere geval van een homogene vergelijking met 


constante coëfficienten geldt 


(u5,17) Vs = Av: 


1+1 1 
zodat 


(u5.18) vs = Alvo: 


i 
Uit (45.8) volgt: 


(H5.,19) De oplossingen van (45.17) zijn dan en slechts dan begrensd 


als alle wortels binnen of op de eenheidscirkel liggen en boven- 
dien de wortels op de eenheidscirkel enkelvoudig zijn. 


(45.20) Opgave. Ga na dat gat (41.5) volgt dat voor een willekeurige 


matrix A geldt dat AÏj danvbegrensd is als alle eigenwaarden binnen of 
op de eenheidscirkel liggen en bovendien de eigenwaarden die op de 
eenheidscirkel liggen uitsluitend 1 x 1 Jordankastjes hebben. 


(U5.21) Opgave. Dat in (45,17) alle wortels op de eenheidscirkel 


enkelvoudig moeten zijn terwijl dit in (45.20) niet hoeft zit hem 
hierin dat in A volgens (45 .M) bij elke eigenwaarde, of deze nu 
enkelvoudig is of niet, slechts één eigenvector behoort, en dus 
ook slechts één Jordankast. Ga dit na. 


(U5.22) We beschouwen weer (U5.15), waarbij A‚v; en b; niet noodzakelijk 


de gedaante uit (45,14) hoeven te hebben. 


AnÂn-1"""A4bot Ss 
TRA 
Er geldt v_,4 * br tAnPn-4 br & ArAn-1 A AoVo» zodat 
(U5.23) Iv, S max IA Arn-1:°:A;ll Ulvoll + oe 


oSisn 
Ons zal de situatie interesseren dat A; z= A + Qi. 


Hiervoor geldt - an 
(45,24) Stelling. Als lIAÌI < MYén IQ; S K voor i < n dan geldt 
n 


IH CA+Q)Il S M(L+MK) TTE, 
lzi+ 
Bewijs. Â CA+Qj) bestaat uit termen van de gedaante 


A...A Q KL, A Q A ... A waarbij sommige der rijtjes A leeg 
kunnen zijn. Een term met j factoren Q heeft hoogstens j+1 rijtjes 
A en is in norm dus < MÌ*Íxà. Er zijn Gee van die termen. Dus 
n n=i _ . « f 
Um (AQ < zE ADM thd, 
l=i+1 j=0 } 


SU6, Probleemstelling beginwaarde probleem. Afspraken. 


(U6.1) We beschouwen de gewone eerste orde diff. vergelijkingen. 


(46,2) TE = ECE 
voor tE I= l[tostg*Tl en xE IR. 
Zoals bekend zijn differentiaalvergelijkingen van het type 
(46.3) Se = gl tsusulsveruftr 
op de gedaante (46.2) te brengen door nl. te definiëren 
Bie) es CapubjnteesBrar Bbbnnn 
als x(t) = (5oCt),Er CE)... vEt). (Merk op dat dan geldt 


a IGE ZEE, 
Dit geldt dus in het bijzonder ook voor diff. vergelijkingen 


van het type 


(46.4) u) 4 a, cHulPTD +…+ artHu = vl. 


(U6.5) Gewoonlijk heeft een dvgl. vele oplossingen en moet (kan) men 
dus verdere condities stellen aan de oplossing. Bekend is het zgn. 


beginwaarde probleem waarbij men eist dat x(to) = Xos Xo een ge- 


geven vector. 


(U6.,6) Een voldoende voorwaarde opdat het beginwaarde probleem althans 


op een omgeving van to precies één oplossing heeft is dat op een 
gebied Q C I x IR", dat het punt (tosxo) bevat, geldt: 

a) f continu en 

b) f Lipschitz in de tweede: variabele; 


Mb 





d.w.z. dat er een constante L bestaat zodat voor alle t;x en X 
waarvoor (t‚x) en (t‚x) € a: 

(U6.7) IECt,x) =— FCE, S Lilx-xll. 
Deze voorwaarde is niet voldoende om te garanderen dat de bedoelde 
oplossing ook op geheel I gedefinieerd is. Beschouw bv. het be- 
ginwaarde probleem En = x?, x(0) = 1, met als oplossing x(t) = Ed 
Men ziet dat f(t‚x) = x? op een willekeurig groot gebied Lipschitz 
is, maar dat desondanks de oplossing slechts tot t=1 bestaat. 

(6.8) We zullen verder aannemen dat het bewuste beginwaarde probleem 
inderdaad een unieke oplossing x* op I heeft, en dat het gebied 
waarvan in (46,6) sprake was, de grafiek van x* bevat, dus de 
verzameling punten (t‚x*(t)) met t € I. 

(u6.9) Zij verder tntotnh, h = T/N, N een natuurlijk getal. Men 
noemt h de stapgrootte. 


Met x! resp. xn geven we de waarde aan van de exacte resp. een 


n 
benaderde oplossing in het punt t,- 
(66.10) We zullen in het vervolg methoden bezien waarmee benaderde 


oplossingen Xn in de punten tn kunnen worden bepaald. De pro- 
blemen zijn daarbij ruwweg de volgende: 


1. hoe groot is de fout x,=xj afhankelijk van h en nh € I vast. 


Geldt Xn” Xn > 0 voor nh vast (d.w.z. voor een vaste waarde 


van t) en h * 0? (convergentie). 
Hoe stabiel is de methode b.v. t.o.v. afrondfouten of kleine 


De) 


verstoringen in x(to)? (het begrip stabiliteit zal expliciet 
gedefinieerd worden en een belangrijk hulpmiddel zijn ter be- 
antwoording van vraag 1.). 

3. Is de methode wel efficiënt? (Het zal blijken dat niet iedere 
methode met succes op elke differentiaalvergelijking toepasbaar 
is, stijve differentiaalvergelijkingen). 


$47. Enkele eenvoudige oplosmethoden. 


(uU7,1) Een aantal bekende oplosmethoden berusten op de opmerking 
dat voor de ware oplossing geldt 


tn 
(47.2) x? me Ee + ff _ fls,x*(s))ds. 


En-p 
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Laten nu XO Xi sere rk al bepaald zijn. 
Zij L*(s) het Lagrange interpolatiepolynoom van de functie 
s » fls‚x*(s)), op een aantal der punten tos---stn-4- 


Men benadert nu (47.2) door 
tn 
(17.3) Xn * *n-p + ER, L(s)ds 
waarin L een interboRatiepolynoom met dezelfde steunpunten als L* 
is maar met functiewaarden f(t;,xj) i.p.v. f(t oxi) voor L*, 
Enkele eenvoudige gevallen: 
CUT .B) Rn FP Br + hÊCt 4%) (Methode van Euler). 
Op CEn-4st) wordt f(s,‚x*(s)) benaderd door L*(s) = FCE): 
% 3 id 
Op (t‚_ost,] wordt f(s,x*(s)) benaderd door het lineaire interpo- 
latiepolynoom L* op de punten tr t,-j: Wordt vervolgens geïnte- 


greerd van tr-, naar tn: 


(7.6) Bovenstaande formules hebben gemeen dat zij op extrapolatie 
gebaseerd zijn. Men noemt ze expliciete formules. Impliciete 
formules verkrijgt men door in de benadering van f(s,‚x(s)) in 
(47,3) het punt tn (waar f nog onbekend is) mee te nemen. 


Enkele eenvoudige gevallen: 
h : 
(H7,7) Xn * Xre1 * zl f(t, xr) + FCE sn) (trapeziumregel). 
Op LEnatn ll wordt f benaderd door het interpolatiepolynoom in de 


punten tred en t- Dit wordt vervolgens van bat naar t, geïntegreerd. 


pe h ee 
(U7.8) Xn 5 Xp-1 * zo Sf Ctx) + Bf Clr an-1) FCE og s*n-2)l- 
Ga zelf na hoe deze methode verkregen kan worden. 
(u7.,9) Om bij een impliciete formule x, te berekenen moet i.h.a. 


n 
een niet-lineaire vergelijking (of stelsel niet-lineaire vergelij- 


kingen) voor x, worden opgelost. Dit gebeurt veelal met successieve 


n 
substitutie (zie 14.9). Men kan aantonen dat dit proces in vele 
gevallen convergeert met convergentiefactor hL, mits hL < 1. (Merk 


op dat hL < 1 geen voldoende voorwaarde voor convergentie is). 
(47.10) De verkregen benaderingsmethoden hebben alle de vorm 


(47.11) XntagXn-1 Eme A Kak = hl Bf CE) tee etB ECE ener )| . 


(47.12) 


SUB, 


(48.1) 


(8.2) 


(u8,3) 
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Zo'n algorithme noemt men een k-steps methode, algemeen een 
multistepmethode of kortweg een multistep voor k > 2, Voor k=1 
spreekt men van een eenstapsmethode. In feite is (47.11) een 
ish.a. niet-lineaire, recurrente betrekking. 


Er is Één opvallend verschil tussen multistep methoden en 
eenstapsmethoden: 
indien x(to) = xo gegeven is kan men direct een eenstapsmethode 
toepassen, echter niet een multistep. Daarvoor moet men eerst op 
een andere manier (bv. m.b.v. een eenstapsmethode) de startwaarden 
Karper eerie bepalen. Men spreekt van het starten van de multi- 
step en van de startprocedure ter bepaling van de startwaarden. 


We komen hier later op terug. 


Algemene theorie van multistepmethoden. 


Bij het bestuderen van methodes zoals (47.11) interesseert 
ons natuurlijk de approximatiefout ee ook wel globale 
discretiefout genoemd. 
Hierbij is uiteraard van belang 
a) hoe goed (47,11) de differentiaalvergelijking benadert 
b) of (47,11) in geval Bo * 0 wel oplosbaar is naar xn 
e) hoe stabiel (47,11) is t.a.v. kleine verstoringen; deze züllen 
vrijwel zeker ontstaan door fouten in de beginwaarden Kore ee sXj1 


en doordat xj niet precies aan (47,11) voldoet. 


Met a) bedoelen we hoe goed een oplossing van de differen- 
tiaalvergelijking aan (47.11) voldoet. Hierbij tekenenwe aan dat, 
omdat in (47.11) f optreedt vermenigvuldigd met h, men eigenlijk 
moet verwachten dat (47.11) eerst links en rechts gedeeld moet worden 
door h om een approximatie voor de differentiaalvergelijking te 
kunnen zijn. Dit is bv. ook heel duidelijk uit (47.4) waar 


On TXn-1)/h een benadering is voor de afgeleide van x* in tn-1: 


Definitie. Zij 
ze en: & 
Dan heet On de locale discretisatiefout (ook wel truncation error). 
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(UBU) Definitie. De multistep (47.11) heet consistent als rn > 0 


“voor h > 0 uniform in n mits nh <T. 


(u8.5) Definitie. De multistep (47.11) heet consistent en van de 

orde p (kortweg van de orde p) voor een zekere aifferant iaal wrdeljking als 
8 P En £ p het grootste getal 
O(h£) voor h 0 uniform in n, het mensne 


n pm 
mits nh £<T 








(48,6) Voorbeeld van wat er kan gebeuren bij een niet consistente 
methode: diff. verg.: SE z= 1, t € IR, x*(0) = 0 (dus ware oplossing 
# . …_ 
x(t) = t). Multistep: Xn Xn-1 5 [FCE isp) tÊCt, sk) le Voor 
de gegeven dvgl. wordt dit 

Xn Xn-1 * oh 

en er zijn voor deze bijzondere f geen startmoeilijkheden: met xo= 0 
is de oplossing geheel bepaald en wel Xn = >nh, zodat men voor 
willekeurig kleine h met nh = t vindt Xn = St, hetgeen ver afligt 


van de x*(t) = t. De multistep is dan ook niet consistent: hô, = 
3 Xn-1 - Sh ze h-sh=-ah wegens x“(t=t, zodat 


* 
n 


Ön =h, hetgeen ‘niet naar 0 gaat voor h * 0, 

(H8.7) Voorbeeld hoe de orde te bepalen: differentiaalvergelijking: 
ES z flt,x), tE IT, ee = xo met nne oplossing x°(t), t € I. 
Multistepmethode: Xn 3Xn-1 +32 ShECt, ox, ). 

Neem aan f € co, dan x* € HS 

* kad * d * En 2 * zi * B * 1 x* n 2 *t … 
Xn I*n-1 * F*n-2 zE Ctx) * Xn F*n-1 * Fn-2 Xn 7 

: Ra rond t=t,_4) 

= «tr + h xr + O(h? )) aj h * 5 -hX + 
à n- 4 4 Jg n-1 FXn-1 * Nije -1 Xn=1 

h «nn 3 me «nn 2 5 
+ Kes + OCH )) sho? 2 + harig + OChD)) = 
= NE 1, be «tt 21 42 «ur 3 es 3 
z (1 zr) Xn + h(1 3 Deil + h (tE F) Xn + Olh”} = OCH). 


Dus is de multistep consistent en van orde 2 mits f € c), Ga 
zelf na wat de orde is indien f € c{2), fe ct1 


(U8.8) Opmerking. De orde van een multistep hangt kennelijk van f af. 
In de praktijk maken we ons van dit probleem af door f steeds "vol- 
doende vaak continu differentieerbaar" te veronderstellen wanneer 


men zegt dat de orde van een multistep p is. 
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(uU8.9) Opgave. Bepaal de orde van (47.4), (47.5), (47.7) en (47.8). 
(U8,10) Zij o het polynoom 
oz) = zh + az | toe + a 


en o het polynoom 
oz) = Bozk + BZT 4 + 6 
o 1 . . k . 


(U8.11) Stelling. Noodzakelijk en voldoende voor consistentie is: 
pl) = 0, P'C1Ì = G(1). 
Bewijs. Bedenk dat x* € C*(I) en bepaal ê, door rond t = train 
Taylor te ontwikkelen. 


(48.12) Consistentie is niet toereikend voor plezierig gedrag van 
een methode. Beschouw bv. eens de multistep Xn BX S*n-2 = 
z hl B fCt sx) +Bf Ct 4 eXn-1)l- Voor geschikte B, en B, ne deze 
methode zeker consistent. Bekijk echter nu eens de dvgl. Te z= 0, 
tS Mti st. 
Voor de gegeven dvgl. wordt Xn uit Xo = x*(t) = 1 en xj bepaald 


door de recursie 


Xn -6x__4 + oXn-2 z 0, 
Neem aan x4 = 1. Dan geldt blijkbaar 
ln +1. 
Neem nu eens x4 = 1+ €. Dan met (45.9) 
= n= E sn 
hi d ns 


Een zeer klein foutje in de startwaarden geeft dus al zeer snel 

een grote fout in het resultaat. Hetzelfde geldt voor afrondfoutjes 
bij het berekenen van X2,X3s-«« Een dergelijke methode zal men 
instabiel noemen. 

Een minimale eis van stabiliteit is dat een foutje in een 
startwaarde of een afrondfoutje onderweg begaan, aanleiding geeft 
tot een begrensde fout verderop, zelfs als h > 0, als er dus zeer 
veel recursie stappen tussen kunnen liggen. Hierbij is aangenomen 
“dat voor h > 0 ook niet de verstoring voldoende ‘snel naar nul gaat, 
wat voor afrondfouten bv. een redelijke veronderstelling is. 

Dit voorbeeld (denk aan (4U5.19)) motiveert nu: 


(UB8.12) De multistep (47.11) heet stabiel indien de wortels van het 
bijbehorende polynoom p binnen of op de eenheidscirkel liggen en 


bovendien de wortels op de eenheidscirkel enkelvoudig zijn. 


(H8.1U) 


(18,15) 


451 


Tenslotte nog het begrip convergentie. Zoals we al in (48.12) 
aan een eenvoudig geval hebben gezien kan een multistep al dan 
niet convergeren afhankelijk van de gekozen startprocedure. Om 


moeilijkheden hiermee te vermijden definiëren we: 


Definitie. De multistep (47.11) heet convergent indien voor elk 


beginwaarde probleem waarvoor (46.8) geldt en iedere startproce- 


dure, die voldoet aan 


lim Xs 5 x(t) ke Tudors erl 
h>0 
de oplossing {x‚} voldoet aan 
Aim ex = x*Ct) 
hn>0 
to+nh=tEI 


Stelling. Voor convergentie is nodig stabiliteit en consistentie. 


Bewijs. 

a) We bewijzen eerst convergentie >” stabiliteit. Zij de multistep 
niet stabiel. Beschouw de diff.vgl. x' = 0, xj * 0» zodat x*(t) = 0. 
Er geldt nu x,ta,Xp.gt: taj XnekK 5 0. Nu heeft de recurrente 
betrekking u‚tagur.gte:-+ Urok 5 0 zeker hetzij een oplossing 
waarvoor lu, |= n (nl. als er een meervoudige wortel op de eenheids- 
cirkel ligt) hetzij een oplossing lul = u" met u > 1. In beide 
gevallen is ook Rs u, /h een oplossing van de recursie en er geldt 
Kres Xed > 0 voor h » 0 en Xn > @ voor h + 0 en nh +» t. 


Derhalve geen convergentie. 

€) ltr Aat pi)=0: Beaehruur xi-0, A3 =1, Zvdat- 
xli=1. de Weurau Áo) ma Xn Al Amert think =O. 
Me Hounds Ho Am =d vor hek DL 5 U 


toegdaAen Wantpreedunrt. _ de Ul aub Lgf Mu tad” 
haf, Ko a Eel a Gar 2 ae E/n, 
Hek Arne. Aa. laat Ky PAAr lt 1, 


Jelle ve dl Verre gef 13. —Ö. 


152 


e) Auw dot p'(1) -0(!), rul de welbedekag Aat zee 
pl)=0 Bee 4) oe PI) golzia). Bogehmus HA, 
„0 zodat alet . A gld he Kn tbAnjt + 
oink = Lo). Bos ehotus- Mu de We ke nb e(1)/pl1) 

vere gh, A Be Mraardi Ho, Aj, VIANEN Ale ole 
ges von Llarfrraards.. Lenh glAf Km +k Ant - Fl 


Â TE | tte) zkronlle) + telle | en hb o(1) 
=p(t)j=o0 —P'0) 

zodet au. de Peewee ò vrtdaa Kuowent eb tm, 

me 12,2 … Dan gld An= EOO/P(I) ven AA oet 

haar t bow Nn oo ‚ Ôeo el) -p0)- 
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$49, Convergentie. 


In deze paragraaf zullen we de fout Xn at nader bezien. We 
zullen voldoende voorwaarden geven, waaronder de fout naar nul 
convergeert voor h * 0. 

We veronderstellen dat de verkregen numerieke oplossing (zo 


deze bestaat, wat we nog moeten bewijzen) voldoet aan: 


(HI,1) Xt AyXrg Peret AKK 7 hl BOfCE sxndte.-tB FCE jen tE, 
waarbij in en de afrondfouten kunnen worden opgenomen en de 
fouten die het gevolg zijn van het niet exact oplossen van 
de vergelijking bedoeld in (47.9). 

De startwaarden KprXgee ee rjg Mofen teGeVe KOX red 
lichtelijk verstoord zijn, dus xj 5 Xi + e, (i = O,1,...sk=1). 
Daarnaast weten we dat de exacte oplossing voldoet aan 


X* HAKS „te. eta, X* 
n 1*n-1 had n=k 


ë At Dan voldoet te} aan de recursie 


eh BGECEr xd teeetB FCE hek) *hôr- 


Z1J en 


ze * En * 
(49,2) eta, egte -td On-k 5 h[ BjfCtnsxnten) BofCtnsxpdteer eet 


, « dn » ee, 
5 Bf CEnejrnek* nek? Bf CEn-k*nek) En hé, 
met de startwaarden eg» esters 


We beschouwen het eenvoudige geval x ER. 


W9.3) Stelling. Zj de methode conscstint . 
0 1 eo \ 


(a) Zij A= | 
Le, 


\za, . ee „al 
(ii) Zij HAI & M onafhankelijk van Nn, Met. 


(iii) Zij h ed (L is de constante uit (46.7). 


28 L 





(iv) Laat een rij Ersree sE, Begeven zijn zodat voor alle 
n < N geldt dat CEsXn + e) € U voor alle e met 


n 
(wv) lel € u cme MYT max le.l+ Elies! + Tô ) waarin 
Osiek, jd ik Hi max 
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wi) ys EK ig a,;l +Je,l) 


2 L 
(vii) c 5, Cha, 1+ ht B; |L) 


(viii) en rn: [6 
Dan bestaat er een rij Keer sXN die voldoet aan wo 
Voorts, als Kore re kid voldoen aan (49,1) voor 
n=k,...,i=1/ dan is (47,11) voor n=i oplosbaar naar Xn 
en wel met de methode van successieve substitutie, die 
zeker convergeert binnen. de grootste "bol" B (dit is in ons 
ééndimensionale geval een segment) met Xn als middelpunt, 
zodat (tu) € 2 voor alle u C B, 

Voor de genoemde rij {x,} geldt 


s me MyT n 
CX) Vxn xrl < 2Me (max le;h 5 le; l+T Ömax?” 
kSnsSN. 
Bewijs. 


We nemen voorlopig de existentie van de rij {x;} aan. Dan volgt 
vit (49.2), als we nog schrijven f(t, oxmtenp-f (t‚sxn)= 


3 Ld 
1, e, ‚ waarin 1, wel van e, en Xn afhangt, maar voldoet 
aan (1 IL, 
(ug, 4) (A-hB len +e--tlar hi Iny) nek” Et Pé * 0 
Zij nu 
en 9 
pek 0 
(44.5) VARA { s b= s 
oe id Ee, 
ë 
n=1 
: 1-hfl, 


(49.6) A= | G Sad 


—a. +hB, 1 - 
Kk RBidnek ar thi ln (wel. 4544)). 
1-hBol, l- hêo 1, + 


se 
[ met “ie Xt EC voor Üzog...k-4 


155 


Dan geldt 
CD) äner * AnZn * Pm 
= b, + AnPn-it ee etAn Ane Are 
GE is de vector der fouten in de beginwaarden) zodat 
n 
49.) Hz, IS max WA ...AL ( Hz,h + EZ Hb.l). 
n+1 keepen n D k en 1 


Schrijven we nu Ans A + Cn dan geldt met de in (vi) genoemde 
Y: 

(as g) HCl, Shy. 
Derhalve geldt, wegens (45,24) 


f n Myhn 
(4g, 1) J VAgsss Al, S CSM) M SM e 
Dus E 
Myhn 5 ‘ 
(43.11) Iz sile SMe C hz, + Pe bil): 
€; 
Nu geldt Ip! nn waarvan de noemer wegens (ii) 
1-hB lj 








minstens } is zodat 


n 
(4q 412) fel & UZ aas S 2 M eliYan, bz hoo + ie 
n 
p 
k 


Wegens hn S< hN = T volgt hieruit ( ix) 

Blijft nog te bewijzen onze in het begin gemaakte aanname 

dat de rij (x;} bestaat en dat CE sx) ER CnsOsletyaee ds 

Dan is (1x) al vast Waar met m i.p.v. m mn % mh ühv | 
We schrijven (49.2) als e, = $(e)- B 

De "bol" B voor X gaat nu over in een bol B voor e‚ 

met O als middelpunt en als straal minstens het rechterlid 


van ( v). & is dus gedefinieerd op de bol B, 
We gebruiken nu de volgende vaste puntstelling,ean 


omformulerine van NÂIQ@6.b6)voor het geval kz}, 
(ug.15) Zij & gedefinieerd op een bol om 0 met straal r, 
zij Id) Pv) $ 3 Ixeyl voor alle x,y € B, en zij 
I6(OL & Ir. 
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Dan heeft de vergelijking $(x) = x een oplossing binnen B, en de 


methode van successieve substitutie convergeert voor elk startpunt 


binnen B naar deze wortel. 


Aan de eis | ò(x)-$(y)| < J|x-yl is in ons geval zeker voldaan 
wegens (iii). 


Voorts is 


(H9, 1u) HLO) = Cray +thByljg)en-it: tra thBj In en-ktEnhên 
en dus 
(49.15) [&o)| S e.maxle;|te,-hö, (ce als in (vii)) 
i<n 
MyT les 14E le, | 
< c.2Me (max e, |+ e: [|+(n-1)hé )te_=hô 
0Sisk-1 Ì iek Î max n n 
ia ar OS ec. 2Mel YT (max le; |+ 5 le; | +nhômax) 
geeen 0Sik-1 isk 
ez (Consistehià) E 
&< halve straal van B. 
Hiermee is alles bewezen. 
(U9.16) Het algemeen geval x € R" kan geheel analoog behandeld worden. 
We doen dit hier niet. 
Opmerkingen 
43.3) 
(49.17) Een multistep die aan de eisen van stellingY“voldoet is blijkbaar 
convergent (zie def. (48.14)). 
Aan de eisen van de stelling is voor voldoend kleine waarden van 
h zeker voldaan als a) AFI < M onafhankelijk van n, b) er zijn 
geen afrondfouten,; en c) Brest > 0 voor h » 0, maar dat zijn juist 
de eisen van stabiliteit en consistentie. Derhalve ‚ samen met 
(18,15) de belangrijke eigenschap: 
(49.18) Stelling. Stabiliteit en consistentie zijn noodzakelijk en voldoende 


voor convergentie van een multistep. 


(u9. 19) 


(u3.,20) 


(49.24) 
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2 behoeft niet groter te zijn dan het gebied beschreven in (iv), 
zodat in feite slechts eisen aan f gesteld worden op een 


onmiddelijke omgeving van x*. 


Voor begrensde T is het effect van een enkele fout e, begrensd, 
hoe groot het aantal stappen (d.w.z. hoe klein h) ook is. 


Ì wd startwaarden. 
Hetzelfde is waar voor de fouten e; in de 


Neem aan dat voor h > 0 de discretisatiefout O(hP) is. Dan 

geldt m.b.v. van stelling (43.3) in afwezigheid van afrondfouten 
(dus alle €; =0) dat de approximatiefout o(nP) is. 

Men interpretere dit niet verkeerd. Zij x Cen) resp. Xn/2 rn 
een benadering voor xXXt) met de stapgrootten h en h/2. Stel- 
ling (49.3) zegt nu dat de afschatting voor DE xp Ct) | 

een factor 2P kleiner is dan de afschatting voor CE) xn CED 


(onder de aanname e; = oChP) en Ee, 0). Dit betekent niet: 
e Pp 6 
[xCt) Xn/2 Eon) |v2 |t) xnCt)l- 


Dat in de praktijk dit laatste toch redelijk blijkt te kloppen 


kan wel aannemelijk gemaakt worden. Men moet dan bewijzen dat 


Xejn “HUE = nP ett) + ochP*Í), 


waarin e(t) niet van h afhangt. Voor eenstapsmethoden kan dit 
in het algemeen wel aangetoond worden. Voor multisteppen ligt 


het veel moeilijker. We bezien dit in €50. 


(50.1) 


(50,2) 


(50.3) 


A58 


De nu klassieke foutschattingen en convergentiebewijzen zijn 
voor het eerst gegeven door Dahlquist (1956, zie ook zijn dis- 
sertatie (1959)). De theorie van Dahlquist wordt zeer breed- 
voerig uiteengezet in de beide boeken van Henrici (zie litera- 
tuurlijst). Ansorge - Hass breidt de theorie verder uit naar 


partiële en speciaal hyperbolische differentiaalvergelijkingen. 


Nadere Foutbeschouwingen, asymptotisch gedrag van de fout. 


Doel van deze paragraaf is een beter inzicht te verkrijgen in 
het mechanisme van de foutvoortplanting bij multistepmethoden. 
In het bijzonder zal in vele gevallen een verscherping van de 
foutschattingen in (49.3) mogelijk zijn. 

Dat dit wenselijk is, moge blijken uit de 

Opgave. De differentiaalvergelijking mX = = Xx wordt opgelost 
m.b.v. de multistep (47.8). Neem aan: x(0) = 1, eg Fe, = 0. 


Jh 


e= 0 alle i. Laat zien dat (49.3) een relatieve fout van 


hoogstens 1% in x(12) garandeert indien h < 10”®, 


In de praktijk zal voor het probleem uit (50.2) een stapgrootte 
h v ta al voldoen. De schattig (4u9,3) eist dus dat de hoe- 
veelheid werk ongeveer 10.000 keer zo groot is als in werke= 
lijkheid nodig is voor het geval uit (50,2). Dit alleen al 


motiveert zeker een nadere foutbeschouwing. M zelden 


daarom Mu Daehl6 ht La bovrtunstua (zaate wo (49.3) 
Maan Us Lbh Pekal Dt Aa fout Ú tes boor 
Ielhotnrt luis hraandter U Á. 


4 


(04). Ze de mealhiakp au orde prin do 450) +0”) 
Me nomen aan dat En (ad (991) Kla boord) Loud, at. 
AO: Er ÔAP*) ME noms veh aan Aat de 
Hanh feude E bolo aa U e= OCAPTL) rekt 

Da. Hert (ega) les dat en = CAL) bro 
Magen we wo (49.2) Alles atrer)- Flbc AL = 
(ÉL ‚) en mat Be (ndrdtee) Aramgpe Let 


(El) es + ÔK ec. h dan. (0 lee) 
Á geldt dus 


Kr 
Escers Ape do ens 


edt Den DAP) az opgponte LE, 

he dk ba AAO nd Ie 

beper). 

Degene es HAL) Moer he te Ôm e AT 
dat dus be behuhe. dOr Uy bengrertglad ls Are 
oo mey, beloken dat 12 PPD ALD LOA) van vak 
4. dat Ape lt), dao hoige he Gat Aonsl 

(50,6) 2 Te Vk NCR: 2 ALs 004) 0 
Bora nalahä deet fork draken Alu Alb ofLedares pat” de 

steady Cevehruurde madhidep Ie de Vngelhoig 


av / eee G „+ 
Gor) Ea opl do Bet be 


(so 5) le A ÓP") zo 
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par valt (sv,6) mat (491) tamter A0 We 
- ÔCA?) noms (NB Me Bedoel Aus Ha 
rn) d Ój (sor) fatende) 
Mo pater te berwoor (49.3) De | Me Aaggen 
da, uit (Hp 3) ex): ne: 
(nl) | ze „ale ame UT max le is, atl+2, DD TA Le, ) 


boachaurn me (TM) het bnaarnerdt ZE ö d 

gat ÈF OA) mor ceker, Orijf Kracht de 

Ahura tes lat. boas De — OU) boor cs ker. Áohabre 
(4) Zr OA). 


orla df Bonte) Alle, we 
(2) vn Kn at - ALU) + OCP) de CU)=0. KR 


Na Hoek Ê* pel de Armer fat Hi edn 
b O/T) de foutemalide (Oror cralat) hea de 
Gon) Aga Kohn onden de aavanaue dat Aa OL) 
mn t (nlo) verwoardoeard Mag wed, [503) Ha. 
CA 5 Wit Murad a Anr alpa ien 
bealoelde ves trdendland Kle do; be de Jrahlyk 
ler zeil Arlrrnnlituk Verbouwen Ab U 


NVrachafft 
Arfa, dg buik Obs waarde Áo, blj Make 
zer / 5 Ag ul VARTA En on) el AV 40) 
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Lorin Alf (A) de met dlp ho, berkraogn. bovadorde 


Plradug von, ira KLE) bee UA) 
(0) Gmerkig — B of) hom EN in ale 

boetaa. drarnm Kum we Hesond (507) allle madr 

gerda ab «* wr Épri)- A Afgelesoa hup De 

bererr ech mel mat meiden (e. TU/ a voldoende 

da nt OD afs antdetaman Magen we dan DE) 

de We Ld HA Jecht j haar INL Mph ot 

dan schriza, TAUEL , Zerbo mlv vat anrboog aa (AÍ). 
os) Vomerkeeg … Úl (507) widt Mor (vet da Ag. meerde 

Ia Maril Me Cella) 
ler) WG f 

(1) 4 
de deffrtud aalndljf rl lb) Ha Wes bedr 
Ar v- Uli) Krt gl A De (Je CIS), ba LAUS ob 
é 

folie Li L ECO do ploaag ha de 

Lrmopant deffrraadbaahren plaag x= CIK ben 

où zò Abfand TA waarde 1 hef don nu te 

(504) de v/ É vewanger doer td, Minten. don 

La GP) ), geduld door ol), eik bofdaa 


É 
é SL afds ie 
EG d id np) AT 


dus NRU vrortsopbear; 


Go.) a (50,3) 





(5u. 1) 


C5H.2) 


(54,3) 


(5 .U) 


(54.5) 


(5U.6) 


em, 


ë 


ee 


Er 
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Beginwaarde problemen. 


Eén van de meest eenvoudige beginwaardeproblemen, 


geassocieerd met een patiëêle differentiaalvergelijking 


is wel: 


t Lex 
ut) sf} 


Gevraagd een oplossing van (54.1) voor t > Os 
De differentiaalvergelijking uit (54.1) staat bekend als 
de warmte- of diffusievergelijking en is van het parabo- 


lische type. 


We veronderstellen f periodiek, periode 1 en f 5 ge 


De Fourier-reeks van f convergeert dan zeker uniform. 


Het ligt voor de hand (54.1) m.b.v. differenties te 


benaderen op een rooster, maaswijdte At in de t-richtin,. 


Ax in de x-richting. Dan is 
B ee ME A U 


J datt J he 
At (Ax)? 


een benadering voor de differentiaalvergelij 


ct 
D 
Uk 
Eh ie 
jan 
lan) 
@) 


(54.1). Hierin is U de waarde van de roos 
in (jAx, nAt). (54b.4) noemt men wel een differentie- 
schema. Een ander zeer bekend schema is dat van Crank 
en Nicholson: 


je + ie Rn me 
n n+t1loyeti + UÈ & Ui -zub +uÏ 


At (Ax)? EAmJS 


Op de gebruikelijke wijze wordt consistentie gedefinieer 


Zo b.v. voor (54.4): Men substituere u(jAx,nÂt) voor Ús 


met u de oplossing van (54,1). 


dd. zr De 





LO 


(54.7). 


(54.8) 


(5uU.9) 


(54.10) 
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Definitie De discretisatiefout van (54.4) in 


(jAx‚nAt) is: 


Made ui -2uk + ui 
DE ae 1+4 ) Jed 
J At (Ax)? 
(5.4) heet consistent indien lim sf = o, nAt en 
Ax0 
jAx vast. At>o 


Nu is het duidelijk dat voor u € C* geldt 

85 = O(At) +0({Ax)2). 
Geheel analoog wordt de discretisatiefout voor het 
Crank-Nicholson-schema gedefinieerd. Evenzo consis :… 


Opgave. 
Ca na dat de diser.fout van Crank-Nicholson 


OC (AH) + OC (Ax?) is. 


Aanw: ontwikkel in Taylor rond ( jAXx , (n+3) BER, 


Zoals altijd is het probleem convergentie aan te ton: 
d.w.z. te bewijzen 
lim us z= ulx,t) nAtzt , jAx=x 


Ax>o 
At>o 


(waarbij i.h.a. tussen Ax en At een zeker verband 
bestaan). 

Evenals bij gewone eerste orde differentiaalverge” 
speelt het begrip stabiliteit een beslissende rol. 


Nemen we voor het vervolg aan Ax = 5 ‚ N natuurlij” 


Door te eisen dat de roosterfunctie eveneens perpiatie” 


is met periode 1, dus Ui in=t5 voor alle js kunner we 


alle differentievergelijkingen (54.4) samenvatten „3 
ARI n 


. n . Ee . 
waarin U" de vector met N componenten is, met als 


ad En 


e . n . « 5 A 
coördinaat Us 4: C is een matrix, die we nu niet 


zullen bezien. 


(54.11) 


(54.12) 


(54,13) 


(SU .14) 


(54.15) 


(54.16) 


NA-165. 


Opgave Ga na dat het Crank-Nicholson schema te seni 
is als 
+ 
Aut*l su? p 
n+1 n ES ER, 
(en dus ook weer als U == ú, indien A nohsijul: 


met een van. (54.10) verschillende matrix C). 


Definitie (54.4) heet stabiel indien heFiezie (5u.ic” 


uniform begrensd is voor alle n en At met o “ nAt # 


T > o. 


za 


Algemeen heet een schema stabiel indien het op de 


gedaante (54.10) gebracht kan worden met hef unite 


begrensd zoals in (54,12). 


We onderzoeken eerst de stabiliteit; daarna zulieu 


aantonen dat voor die combinaties van At en Ax waar: ” 


(54,4) 


stabiel is ook convergentie optreedt (en cv 


voor Crank-Nicholson). 


We maken nu gebruik van de periodiciteit. 
Zoals bekend kan de exacte oplossing geschreven uus, 


als een Fourier-reeks. Dit doen we ook voor de on! 


op het rooster. We stellen 
N-1 En 
lóieh s 5 ar Go) aikTj/N 
J _kzo 

Er geldt: 

ee ik j/N =imrj/N _, O Kk #m 

: e EE: ={ à 
j= N k =m. 
j=o 


km == OslsèsgeessN=1 


waaruit volgt: 


N=-1 NEN 
5 De er Her/N, 
in 


uwék) = 
o J 


zi 


We kunnen dit ook als volgt interpreteren: def inie 


lineaire ruimte van de functies gedefinieerd op de 


J/N, 


jzo,l,s...,N-1 het unitaire inprodukt: 


(54.17) 


(54,18) 


(54.19) 


(54,20) 


(54,21) 
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N-1 


24 A 


(EE) = FCGÍN) ECH/ND. 


zin 


j=0 


Den de uiikh es B ; sj terwijl (54.15) aantoont 


LT 217 N- LT 
die oiTx e2i X, EEE ml 1)irx 


basis t.a.v. dit inprodukt vormen. Dus 


een orthogonale 


n Bea 
k 
en (UP) es UP s E (UP, eIKTxj2, NEN % (GG)? 
k 


Deze laatste relatie staat bekend onder de naam Parceval- 


relatie. 


Substitueren we (54.14) in (Bu.k) dan vinden we voor û'(k) 
de recursie: 


A 


uto) = elk, aP Go 
met als amplificatiefactor G 


GCK,AE) = 1e u. AE)z sin? ar 


At À 
“(Ax 2N° 
Stelling (54.10) is stabiel indien le? (k,At)! 
uniform begrensd is voor alle n en At, o <nAt < T en alle 


Ni 1 n 
Re. HU te KU H. 


Re _ as 
B Ons te wr nn ee arr te sE 
j dk 
UL, UD =N.r [eFk,At) a° (12 
j J k 
&M.N.EB Hi Cleut * S MELLE U jh. 


k 
indien M de uniforme bovengrens voor |G (k,At)lis. 


Dus: 


n 
Eer < M waaruit volgt ICI, SM. 
0 


(5U.22) 


(54.23) 
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Opgave Ga na dat de voorwaarde van (54.21) ook nodig EN 
Voor (54.20) betekent de voorwaarde van (54,21) kennelijk 


At < 1 
(Bah 4 
Onder deze voorwaarde (dus stabiliteit) volgt eenvoudig de 
eerder aangekondigde convergentie van (54.4). 
Zin mm. el de vector (N coördinaten): 
n h_f 


e = U =u 


de 


: 5 add n 
waarin de j coördinaat van u gegeven wordt door 


utd en nAt). ef stelt dus de fout per horizontale 
rij voor. Dan volgt uit (54,10) en (54.6): 


n+1 n 
e 


stel + At‚sf 
waarin 81 de vector van de discr. fouten voorstelt. 


Dus als ICI, < M: 


en+1 = cr+1 e + At{ cls: + hie PRE 5] 


n+ 1 


le 2 S Mell, + At.M.n. max H8I, 


< Mlle® ll, + M.T. max Il8 Il» 


waaruit de convergentie volgt. 


We hebben nu convergentie, maar niet in de sup=-norm. 
Ô : de n Dn 2 
Maken we gebruik van de periodiciteit en f 5 c{ ) 


dan 
impliceert convergentie in de gebruikte norm ook uniforme 
convergentie, d.w.z. convergentie in de sup-norm (c9- norm). 


Maar we kunnen niet zeggen dat Let „ 5 Olmax hs). 
Ä k osisn cl 
Dan zouden we moeten weten of I|C | ie, uniform in n 


begrensd is 


(5u .24) 


(54.25) 


(54.26) 


(SU .27) 
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Op het Crank-Nicholson schema kunnen we onze theorie 
ook toepassen. Zonder eerst na te gaan of A 

(zie opgave (54.11) non-singulier is, substitueren we 

weer (5h.14) nu in (54.5). 

Dan vinden we voor de û'(k) de recursie: 

_ 1-2 Te: sin? sr 


kr 
1+2 )2 sin DN 


Le) ak) 


(Ax 


zodat de amplificatiefactor aan de eis van €54.21) voldoet 
voor alle Ae? ‚ Dit impliceert dat het Crank-Nicholson- 


schema een unieke oplossing bezit en dus dat A nonsingulier 


is voor alle Ai)z . en ook Art B, uniform in n 


begrensd is voor alle n. (o <nAt<rT). 


De beschreven methode is natuurlijk alleen toepasbaar 
bij een differentiaalvergelijking met constante coëfficiënten. 
Maar dan loont het ook wel de moeite. Bezie b.v. de warmte- 


vergelijking in twee dimensies: 


ie Aex ij Uyy 


waarbij als startwaarde een functie f(x,y) gegeven is, 
periodiek in x- en y-richting met periode 1. Het is eenvoudig 
een differentieschema analoog aan (54.4) op te stellen en dit 
kan ook nog in de vorm (54.10) geschreven worden. Analoog 
Crank-Nicholson. 

De optredende matrices C zijn echter bijzonder onplezierig 


terwijl de methode m.b.v. Fourierreeksen eenvoudig werkt: 


Zij nml. 
oat _ n n tn n n tart n 
dakar etiket lier Pier Piet Dies 
At (Ax)? (Ax)? 
n 


het analogon van (54.4). U. k is de waarde van de rooster- 


b] 
funetie in (jAx, kAx, nAt), Ax=d d 


Stel nu UL z= % â(p,a) e 


(54.28) 


(54.29) 


(54.30) 
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ipTi/N 8 iqTk/N 


k k2 
SET Bad 
We vinden voor â(p‚a) de recursie: 


ant1 n … uât oat BL At ral 8E ss 
A Pes Cl Wipe OD mT bere od Â pad). 
en geheel analoog aan het voorgaande is er stabiliteit en 
convergentie indien 
u ät tat BE … At ent AL A 
eb Biz Sin? Sp bram Sin? pn) 
! 
uniform in n begrensd voor alle p en q, dus (1x) 2 hd 
(Merk op dat de verhouding tussen At en Ax hier veel 


ongunstiger is dan voor het schema (54.u)). 


Men kan ook de methode toepassen op niet parabolische 
differentiaalvergelijkingen, b.v. de (hyperbolische) 


golfvergelijking ( in Één dimensie): 


Utr 7 Ux” 


met voor tzo gegeven u(x,o) en u„(x,0) die beide weer 
periodiek, periode 1 verondersteld worden. 
We voeren weer roosterfuncties in op een orthogonaal rooster 


met maaswijdte Ax in de x-richting (Ax = En N natuurlijk getal) 


en At in de t-richting. Als differentieschema kiezen we: 


n 


n+1 
E j+1 


n n=-1 n 
U 2U. + U. U. 
Ki EREN et ROE oct MIND | 


CAE (Ax)? 


au HU 
ki 


In tegenstelling tot de warmtevergelijking is de golfvergelijking 
een differentiaalvergelijking met een tweede orde afgeleide 

in de t-pichting. Dit komt in de benaderende differentie- 
vergelijkingen (54.30) tot uiting hierin (naast het feit dat 
(5SH.30) een multistepmethode is) dat (54.30) niet begrensde 
oplossinsen bezit voor nAt€ [o,T] met wel begrensde start- 


waarden kies maar Ui zo en ui = 1 voor alle j en taat 


At > o gaan. We zullen dus geen stabiliteit in de zin van 


(54.10) mogen verwachten. 
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We omzeilen het probleem door (54.29) terug te brengen tot een 
stelsel eerste orde partiële differentiaalvergelijkingen en 
(54.30) tot een differentieschema voor dit stelsel. 

Van dit equivalente schema zullen we stabiliteit en dus 
convergentie aantonen (en dus geen stabiliteit van (54,30)), 
Tenslotte tonen we wel convergentie van (54.30) aan, uitgaande 
van de dan al bewezen convergentie voor het differentieschema 
dat de eerste orde differentiaalvergelijkingen benadert. 


(54,29) is equivalent met 


V 5 We 
(5u.31) 
Ve p he” 
door te stellen w = U, en v =z Ur 
Het differentieschema (54.30) is equivalent met 
n+1 n n n 
u.32 Vv. - Vi s - W. 
(54 .32a) 5 d - Ws PE Ta 
t X 


(onder WE, PE g 


door te SKL 


ee EN: SN 


j At 
us gl 
we = ee ON 
jj Ax 


n 
Zij (-Ímjde gepartitioneerde vector met als zie coördinaat 
W 


de coördinaat van WÌ : w@ 


van VÒ ; Vv DE 
Ei 


B, : 
5-1 en als j 


Dan is er een matrix C zo dat 


n+i n 
(54,33) dels. | gE 3 


weti we 


NAT 7 


We tonen nu eerst aan dat onder zekere voorwaarden (54.32) 


een stabiel schema is. 


We schrijven de oplossing van (54.32) als: 


(54.34a) vi = Yr 


n E: 
(54 .34b) We, s E 


Ne) 8 ik j/N 
& Aa) a IKT (+5) /N 


geheel analoog aan (54.14). Nu is (54.34) oplossing 


indien voldaan is aan de recurrente betrekking. 


À 


ant 


Î E í \ 
iv do, edo | 
le = G(k,At) | 
(58.35) | ntt! Goo | 
waarin G, de amplificatiematrix egeven wordt door 
t id: 
Ea ia : 
=| „8E sin 
(54.37) Stelling (54.21) is in iets gewijzigde vorm toepasbaar; 


het verschil is dat we nu vectoren i.p.v. scalairen hebben. 


Omdat het 


bewijs van (54.21) uitsluitend op de Parseval-relatie 


berust ‚ bezien we deze relatie nog eens voor ons geval 


apart. Daartoe herschrijven we (54,34) als 


waarin ( 


B : 
Ve - ae ue, ekx) 0 ik j/N 
J e 
k=zo 
N-t 8 é n 
Wie zz yr, IRT EKAIN 
J*5 kzo 


‚ Jen < > inprodukten zijn, 


( _„ ) gedefinieerd als (54.17) en < > als 


ze 


Ne tb 5 
TE = 5 EG) gr). 


j=o 





(54,38) 


(54.39) 


(54.40) 
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Dan luidt de Parseval-relatie: 


Dee] 


E dvtaol? + Leao DD). 


ve 2 rij ° N22 … 
HC z VI + IW IS == NC 
2 2 2 k=8 


W 
Het analogon van (54.21) blijft gelden met 
IGP (K,At)ls uniform begrensd (i.p.v. Le txK,At) 1 ) 


en met als matrix-norm van C de geassocieerde van de 2-norm 
(in de ih pe ms, 


Om de stabiliteit te onderzoeken passen we (45.18) op 
G(k,At) toe. Nu worden de eigenwaarden van G gegeven door 


> 2-a? * va“-ig? 


À 
1.2 
2 
Voor azo geldt G = I (2x2 eenheidsmatrix). 


Dus is (54.32) stabiel voor a € (-2,+2), d.i. 


at <1. Voor Ee = 1 is de Jordan-normaalvorm van G niet voor 


alle k een diagonaalmatrix, dus geen stabiliteit voor ie z 1. 


Analoog aan wat voor de warmtevergelijking is afgeleid vinden we: 


BE <1 
X 
n o sP 
v Sy Vv 
IN =d < M. WC-sorlle +M.T. max Wz) 
w p Oer 
er 
waarin (--Z--) de vector van de fouten op de n-de horizontale 
ef n 
w Ö 
rij voorstelt; amd die der discretisatiefouten. 


wW 


Dit toont nog eens de convergentie van het schema (54.32) aan; 
uitgaande van de stabiliteit van dit schema. Als de exacte 
oplossing (vw) en zijn partiële afgeleide (vw) eontinu en in 
de x-richting periodiek zijn volgt ook dat de benaderde oplossing 
naar de exacte convergeert in de sup-norm (dus uniform). 

We weten echter niet of de fout als O(At®), voor zekere a > o, 
naar nul gaat. 

We herinneren ons nog: us = vj 


wat een benaderde integratieformule voor U, =V is. 


+ At.vi 
J 


NA=172 


Als nu maar vi convergeert naar v(j/N,nÂt) voor At > o 
z dan is wel duidelijk dat u; naar u(j/N, nAt) convergeert. 


(54 .U1) Een andere situatie waarin (54.29) optreedt is de volgende: 
gegeven u(x;o) en u, (x,0) voor x & [o,L] terwijl een oplossing 


gevraagd wordt op de driehoek met hoekpunten (oso) , (o,L) en 
(EL,p), Ípl < 4 L. Dit geval brengen we tot het vorige terug 
door u(x,o) en u, (x,0) tot een periodieke functie uit te breiden, 
periode > L. Voor dit uitgebreide probleem is convergentie verze= 
kerd. 

Men kan echter eenvoudig aantonen dat op de beschouwde 

driehoek de oplossing onafhankelijk is van de uitbreidingen 

van u(x‚o) en u, xo). 


Dus ook convergentie voor het probleem op de driehoek mits 


At … 
rn 1e 
(54.42) Indien de = 1 kan ook convergentie aangetoond worden. 


Zij Lb=1. Er geldt: 


n n-1 n=? 2 n 
TE Ee - @. T+ (At ê. 
23 j+1 © S5 kar 8 


waarin es de fout in us voorstelt en 55 de discretisatiefout. 


Herhaald uitschrijven en afschatten: 


neel n=2 
lefl <1 Ee Ssgsgp TE et ig LED at? maxlën 1. 
J p=o J=n P pz1 }J P J 
Nu is e; = o voor alle j (afgezien van afrondfouten) 


en zij e, = O(At® ) voor alle j. Dan: 


Pd 


n 1 : n 
Lesl &< n. max le;! + max | $ | 


1 


ee OAT TE) + OC max 18E 1). 


ED 


Als nu u dan O( max KAR z O((AtT)2) 


zodat als a > 1 convergentie verzekerd is. 


(5u. U) 
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Opgave. We berekenen de startwaarden U; als volgt: 
mg 
J 
2At 
(54,30). Ga na dat op deze wijze een U; verkregen wordt die 


u(j/N,At) benadert met een fout O((At)?).(AxzAtz=1/N). 


= u, (3/N,o) + O((At)2). uz” elimineren we m.b.v. 


Sommige auteurs noemen het schema (54.30) stabiel 
voor ke < 1. Deze hebben dan een andere stabiliteitsdefinitie 


(wie b.v. Forsythe en Wasow blz. 32). 

Merk nogmaals op dät we hier alleen stabiliteit voor een 
differentieschema gedefinieerd hebben dat een differentiaal- 
vergelijking met eerste orde afgeleide in de t-richting 
benadert. De gegeven theorie m.b.v. Fourier-reeksen en de 
gegeven stabiliteitsdefinitie is afkomstig van Lax en 


Richtmyer. Zie vooral het boek van Richtmeyer en Morten. 
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| | Ned 
3 me, Beleg aa flh ssnedjeeign 
4 Bedmnielhs 


(11) da a [CE 2) n xe kk. 
d £ 


met a Ck) BES aaven. 
ak cle Nete Gsbastd pa Í op Ad ons 2 
oghried . WL meme deme Bek oe opFoanin wan. (44) ba 
alst aw dla elg is ep TL, ber T € 
We Belen Le scale opdoemen, k vin olen Armor 
be hbs. RT IP We EE gedker werd op Le 
Numerieke nm amen ves ahh p4 nä L. Naaf cle agen | 
Arna Ainden, , maar Hse ons EN JEN Tepper Hen 
deb at lane elen, ep Teka helf bk kermen panter 
bs We Aachen Emek eme es deed mieekp bare, 


2 Nasse vende hein, 


zn) We dedeminen ten lapse es TIN eN zen 
en 2 agdal. Mel Jl vuAP Am deg kaders cle 
AFA ole Oan aman cle ex mate Aap . MNAn Len NSE 


dede. oplanning die Al pak He dem. bn 
dte bek ophaal deelen kenia of> ele et bni 
elat 


bs 
B te el Hike 
Ë | 


jp 


(2.3) 


(2.4) 


(2,8) 
(2.6) 
(2.7) 
(2.8) 


(2.9) 


Nu geldt voor een willekeurige functie g € ct bijv. 


t 2 
n Ee h 
S g(t)dt = hg(t,__4) + 5 g'( En) En-4 < En < B 


t 


(integratie van het Oe orde interpolatie polynoom in tn-1) 


n_l 


Door toepassing op (2.2) komt er 


2 
Ke Ss Ke Ea Hf + h_ x* (EE ) t < E <t 
7 iks -1 n 


n n= n=1 2 n n 


Analoog krijgt men met een kwadratuurformule als erachter aangeduid: 


2 

NR * EN : 0 PL e . . 
KSK th Lan DX (&) (0 gr. interp. pol. in t_) 

* * 1 * * n° „ed + 
ES Bregt zh[ Én + A =  X (CE) (trapezium) 

dn ik 1 * ex Bs e : * 
Ee, 3h[ 3fL_, Era] + hx CE) hesp pel: ) 

n=-1° n-2 

E B m ns MAA : . 

Xh S Xh-g + Mh fBg + 5 X CE) (midpoint) 
5 

* z x* Ar geruft * „ h „«(5) 

Xn Xh-2 * zn ER+BEN tE] 50 * C&) (Simpson) n 


Hierin is Én = ECE» Xn etc. Door in (42.4) t/m (2.9) de sterre- 
tjes weg te laten;en ook de termen met &, worden dan blijkbaar 
benaderde oplossingen van de dvgl. gegenereerd, raar hu de 
waanden Ks, Ymaa, --- Len waarde X be paden. LN 
olen caire laa krwkense zwe. ha neiburdt oadbekilkanike. 
Formules (2.4) en (2.5) duidt men wel aan als Euler resp. Euler 
badande: 

Het valt op dat: 

- de methoden niet alle dezelfde orde van nauwkeurigheid hebben 
Khel C2eB) en (2.9) Xn niet zondermeer opleveren, omdat X 
ook in het rechterlid voorkomt. Deze formules heten daarom 


impliciet, in tegenstelling tot de overige, die expliciet heten. 


Xx) Dit interpolatie polynoom wordt van tn tot t‚ geïntegreerd; 


dat is dus eigenlijk extrapolerend integreren. 


(2.10) 


(21) 


(212) 


(2.13) 
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Bij de impliciete methoden bepaalt men Xn meestal met successieve 


substitutie Ce sok Anr (6.22)) S Bau voor (2. é) : 


Zx an een beginschatting voor Xn (bijv. zn = Xn)" Dan gene- 
reert men xn, en d.m.v. 
n n 


+1 ì ì 
ge 5 Xn-1 ú En(ECt, x0d)) + FCE roger) Voor h klein genoeg 
convergeert dit: bide - hd < nh iss aal (L ele hirchits end 


n } 3 
zodat voor 3h L <q < 1 geldt Enid Sad < A Fn ed als 


q 'n 
Ì > j, dus Cauchy's criterium is vervuld. De convergentie is zelfs 
snel omdat beginschattingen. als ee EH bij kleine stapgrootte 


uiteraard nogal goed zijn. De reden voor impliciete methoden is o.a. 
dat men op deze wijze interpolerend integreert, in tegenstelling 
tot expliciete methoden, waarbij men extrapolerend integreert, wat 


wel onnauwkeuriger zal zijn. 


De vraag of een hogere orde methode ook een nauwkeuriger antwoord 


oplevert, vereist nader onderzoek. (zie hter Voor bgv. Mum. u. i) 
De verkregen benaderingsmethoden hebben alle de vorm 


X + Q +  . « + Xx 


n 1 *n-1 Kk ak 
Zo'n algorithme noemt men een multistep of ook een k-staps methode. 


= h [BÉ tee. + BÉ pd 


Voor k = 1 staat er een onestep. 


ei 1 * * * == * * 

Zij B 5 {xf tak teen + AKA h[Bjfp+--.+B, Én}: Dan 

heet an de locale discretisatiefout (ook wel truncation error). 
2 ' 

Zo is dus bij (2.6) de locale discretisatiefout - on x* (B. 

De reden voor het delen door h bij de definitie van locale dis- 


cretisatiefout is dat (2.12) eigenlijk t.o.v. de differentiaal- 


vergelijking met h vermenigvuldigd is. Dit is bijv. heel duide- 


lijk bij de regel van Euler: x, = xj + hf, oftewel x, - x__4 * h 6 


terwijl een benadering van de dvgl. eigenlijk zou luiden 


(x, 7 Xj-j)/h = f. De locale discretisatiefout geeft dus in 


feite aan hoe goed (2.12) de dvgl. benadert. 


(2.14) Tot dusverre werd het begrip orde in een intuitieve zin gebruikt. 
We geven nu een definitie. 
Als x* voldoende vaak differentieerbaar is heeft _ steeds de 


gedaante enPxr(P+1) + oqnP*1), Men ziet dit in door in (2.13) in 


: REN 
* = « _ Jh De, (jh) ë = « me nt, si 
te vullen Xn-j ° Fn ES Xn + NE ns Xn …… en Ín-3{ Xn-á’ 
1 D " ' 
Xn = Lr Xn + .…. Men zegt dan dat (2.12) van de orde p is, als 


P de agsotsh woar de La Herdt S de aa apnn ggd aat 


(2.15) lik eN vann sds epen les eee se Ln Len 
leo vrm media, Ault orn Kas oes Xe Ü 
bepalen FN Raap eren ale As en LA. 


(2.16) OENE Mebo (2.12) is a Jl dele dla A Om 
ole or ele, P >={ ale addeld 


EE 


k + (k- 4) aes + (k-2) «2 eat 5 = Se Fen +k 


Den Biapnehap ke Amd ziet: óle Te en En 
Lag Schop, 


ä Velada 


(2.1) Ade DL PN mu ela. wench bag le,Ì 
van Len edna en Len Avane ER 4 Nam Lem 
df) med: he olga duzt ne Us oet 
En ha aanrelig gelast aen: de Apidae, 
c U A « 


2) x 


kad 
s Bme Aen Kars 


= PX + gt) 


he me ans Ht iS Aaen indi Ed & 
en Wen he da 
mak: DN | ti (2.4) 


zen ele 


3.2) OE wk GE 
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bad Brake de Ka le 
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Gekat 14 Lp ee Ws Ake hin eld Wi das he biek- 
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AE + ed on + otk Gate 
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er nak go sec ú 
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